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Abstract 

n 

We study the symmetry breaking of acoustic resonances measured by Ellegaard et al. [1] in quartz 
blocks. The observed resonance spectra show a gradual transition from a superposition of two 
uncoupled components, one for each symmetry realization, to a single component well represented 
by a Gaussian orthogonal ensemble (GOE) of random matrices. We discuss the applicability of 
superstatistical random-matrix theory to the final stages of the symmetry breaking transition. A 
comparison is made between different formula of the superstatistics and a pervious work which 
describes the same data by introducing a third GOE component. Our results suggest that the 
inverse-chi-square superstatistics could be used for studying the whole symmetry breaking process. 
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INTRODUCTION 



Few years ago Ellegaard et al. studied the influence of the flip symmetry breaking in 
quartz blocks on their spectral statistics [lj]. Crystalline quartz exhibits a D 3 point-group 
symmetry about the crystal's Z (optical) axis and three two-fold rotation symmetries about 
the three X (piezoelectric) axes; the latter three axes lie in a plane orthogonal to the Z axis 
and sustain angles of 120° with respect to each other. Ellegaard et al.jlj used rectangular 
blocks of dimensions 14 x 25 x 40 mm 3 , cut in such a way that all symmetries are fully 
broken except a two-fold 'flip' symmetry about one of the three X axes. In order to break the 
flip symmetry, Ellegaard et al removed an octant of a sphere with a successfully increasing 
radius from one corner of the rectangular block. The breaking of the flip symmetry causes the 
eigenmodes belonging to the two different representations to interact and mix. As the radius 
is made larger, the nearest-neighbor-spacing (NNS) distribution makes a fast transition from 
a distribution of two non-interacting spectra towards the distribution of one chaotic system 
without symmetries. Full chaos was attained at an octant radius of r = 10 mm, as the cut 
block had the shape of a three -dimensional Sinai Billiard. 

A problem like this is found, e.g., in nuclear physics where isospin symmetry, characteristic 
of the strong interactions, is only approximate due to Coulomb effects j^]. Isospin mixing 
was analyzed by Guhr and Weidenmuller using a random matrix approach J] . They used a 
random matrix model to describe experimental data and to estimate the average symmetry- 
breaking matrix element, i.e., the average Coulomb matrix element. The random matrix 
ensemble consists of matrices H = H + V, where H is block-diagonal, here two blocks 
for the two-fold symmetry, each a member of the GOE; V couples the blocks, breaking the 
symmetry. Level statistics were studied within this model analytically by Leitner jlj. He 
assumed that the probability density of spacing belonging to the same block are unaffected 
by the coupling matrix V. He took into account the interaction between pairs of levels 
belonging to different blocks, yielding linear dependence of the spacing distribution near 
the origin. Taking into account the interaction between more than two levels might lead to 
fractional-level repulsion as recently shown by Backer et al. 6|. While Leitner's formalism is 
a perturbation-theory approach, which is valid only for small V, he has used it as a model for 



tho whole symmetry breaking process in quartz-block data El-Hady et al. 



described 



the transition by introducing a third GOE sequence which increased in the expense of the 
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initial two independent sequences. To fit the data in Ref . [1] , the authors assumed that the 



■ Q, 



initial sequence are pseudointegrable (for definition, see the discussion following Eq. ffl~2]) 
below). This independent- sequence model of symmetry breaking was further elaborated in 
9], where a relation between the fractional level density of each sequence is related to the 
mean squared symmetry-breaking matrix element. De Carvalho et al. 10] take into account 
the interaction of the three level sequences using the perturbation method suggested by 
Leitner js|. 

In the present paper we perform an analysis of data in Ref. [1] using a superstatistical 
generalization of Random Matrix Theory (RMT) {2, [ill 14]. We are motivated by the 
success of the superstatistical approach to describe the final stages of the transition from 
integrability to chaos. A brief account of this approach is given in Section 2. Section 3 
reports the results of the superstatistical analysis of symmetry breaking. We also compare 
the results our work with that of Ref. [8|. A discussion of the results of this contribution is 
given in Section 4. 



SUPERSTATISTICAL RMT 



The concept of superstatistics was introduced by Beck and Cohen [15J to describe de- 
viation of thermodynamic system from equilibrium. This concept has been successfully 
applied to a wide variety of physical problems, including; turbulence Il6[ , plasma physics 
13], cosmic-ray statistics [isj], cancer survival 19], and econophysics |2o| . The application 



of the superstatistics to RMT [21] assumes the spectrum of a mixed system as made up 
of many smaller cells that are temporarily in a chaotic phase. Each cell is large enough 
to obey the statistical requirements of RMT but has a different distribution parameter 77 
associated with it, according to a probability density f{rf). Consequently, the superstatis- 
tical random matrix ensemble that describes the mixed system is a mixture of Gaussian 
ensembles. The joint probability density distribution of the matrix elements is obtained by 
integrating distributions of the form 

1 



P(H) 



z(v) 



exp 



(1) 
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over all positive values with a statistical weight f(r)), which leads to 

where Z(ry) = / exp[—r]Tr(H^H)]dH, 

One of the fundamental assumptions of RMT is that the matrix element joint probability 
density distribution is base independent. This makes the theory suitable for modelling 
quantum chaotic systems. Indeed, the eigenfunctions of a Hamiltonian of a system with a 
chaotic classical limit are unknown in principle. In other words, there is no special basis 
to express the eigenstates of a chaotic system. In integrable systems, on the other hand, 
the eigenbasis of the Hamiltonian is known in principle. In this basis, each eigenfunction 
has just one component that obviously indicates the absence of complexity. In the nearly 
ordered regime, mixing of quantum states belonging to adjacent levels can be ignored and 
the energy levels are uncorrelated. 



Parameter distribution 



The distribution f(rj) is determined by the spatiotemporal dynamics of the entire system 
under consideration. Beck et al. 2^ have argued that typical experimental data are described 
by one of three superstatistical universality classes, namely, x 2 i inverse x 2 or log- normal 
superstatistics. The first universality holds if rj has contribution from v Gaussian random 
variables Xi, ... , X v due to various relevant degrees of freedom in the system. Then a 
positive 77 is obtained by setting 77 = Y^i=\ X? is and / (77/770) is a x 2 distributed with degree 



fWlo) = -^-(f)^ 2 - 1 ^ (3) 
r(z//2) 2r] 

Here, i] is the average value of 77, with 770 = / °° 77/(77/770)^77 

In the second universality class, the same considerations as above can be applied if r]^ 1 , 
rather than 77, is the sum of several squared Gaussian random variables. The resulting 
/ (77/770) is the inverse- x 2 distribution given by 

fiv/r]o) = ^0r /2 ^ v/2 ~ 2e ' uvo/2v ( 4 ) 

In the third universality class, instead of being a sum of many contributions, the random 

V 

variable rj may be generated by multiplicative random processes, i.e. 77 = 'f[Xi. Then 

i=i 



In 7] = Y^=i ^ n Xi is a sum of Gaussian random variables; hence it is Gaussian as well. Thus 
is log-normally distributed, i.e., 

f(v/Vo) = ~f=~ e V , (5) 



'2itut) 

which has an average fiy/^p and variance a 2 = /J, 2 ip(ip — 1), where ip = exp( 



NNS distribution 

It follows from Eq. (j2J) that the statistical measures of the eigenvalues of the superstatisti- 
cal ensemble are obtained as an average of the corresponding ^-dependent ones of standard 
RMT weighted with the parameter distribution / (77/770). In particular, the superstatistical 



NNS distribution is given by 21] as 



p( s ) = f(v/Vo)Pw(v,s)dr], (6) 





where p w (r], s) is the Wigner surmise for the Gaussian orthogonal ensemble with the mean 
spacing depending on the parameter, 

Pw(v, s ) = 77sexp(-^77s 2 ). (7) 

For a x 2 distribution of the superstatistical parameter 77 the resulting NNS distribution is 
given by 

P * 2(M = (1 + 770,7^ ' () 
The parameter 770 is fixed by requiring that the mean level spacing (s) equals unity, yielding 

% = T [r(^-)/r(-)] 2 . (9) 

For an inverse y 2 distribution, one obtains 



Pinv X 2 (^s) = ^^(^ri vs/2y /2 K v/2 (y/ri vs), (10) 



where K m (x) is a modified Bessel function 



and 7/0 again is determined by the requirement 



that the mean level spacing (s) equals unity. Finally, if the parameter has a log- normal 
distribution (5), then the NNS distribution 



— 2 \n 2 (^ne- v2 ' A ) 



PLogNorm(v, S) = exp[ — \dr] (11) 

V Ittu Jo 2. Lv* 



can only be evaluated numerically. 



A justification for the use of the above-mentioned superstatistical generalization of RMT 



in the study of mixed systems, is given in 24( . It is based on the representation of their energy 
spectra in the form of discrete time series in which the level order plays the role of time. 
Reference 24| considers two billiards with mushroom-shaped boundaries as representatives of 



systems with mixed regular-chaotic dynamics and three with the shape of Limagon billiards 



one of them of chaotic and two of mixed dynamics. The time series analysis in Ref. [24 1 
allows to derive a parameter distribution f(rj). The obtained distribution agrees better with 
the inverse x 2 distribution given by Eq. fll]). The inverse x 2 distribution of rj follows when 
the quantity r) is the sum of v inverse-squared Gaussian random variables. In the application 
to RMT, the parameter rj is proportional to the sum of the inverse variances of the matrix 



elements of H 24]. Hence, v refers to the number of (largely) contributing matrix elements. 
In the limit of v — > oo where all the matrix-elements contribute to the distribution, Eq. (j4]) 
yields a delta function for Pi nV x 2 { v i s ) turning superstatistical matrix-element distribution 
([2]) into that of the conventional RMT, Eq. (CE]). If we take this assumption literarily, we 
must restrict v to take positive integer values. As the transition from integrability to chaos 
is known to proceed continuously, we have to relax this condition and allow v to take any 
real value grater than 1. Using the asymptotic expression of the modified Bessel function 



23(|, we easily find the Wigner surmise when v — > oo, as required. The other limit of v — > 1 



yields the semi-Poisson distribution 

p ss (l,s) = 4se- 2s , (12) 
which is known to provide a satisfactory description for the spectra of pseudointegrab 



systems such as planar polygonal billiards, when all their angles are rational with tt 25]. 



The motion of the corresponding classical systems in phase space is not restricted to a torus 



like for integrable systems, but to a surface with a more complicated topology [26(. We 
therefore conclude that the assumption that the inverse square of the variance of matrix 
elements as an inverse x 2 variable allows superstatistical RMT to model the transition out 
of chaos (corresponding to v ^> 1) until the system reaches the state of quasi- integrability. 
It is interesting to note that, if one allows v to take lower values, one finds that the 
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distribution ( !T0|) tends to the Poisson distribution as v — > — 1; 



p S s(-M) = e- s . (13) 

We there conclude that formula (fTUI) can provide a successful model for describing the 
stochastic transition all the way from integrability to chaos passing by the stage of quasi- 



integrability. This has been clearly demonstrated in Ref. 24]. We have no physical ex- 
planation for this success. We regard psg(v,s) in the range of — 1 < v < 1 as a clever 
parametrization of NNS distribution of nearly integrable systems undergoing a transition 
from a Poisson to semi-Poisson statistics. 



DATA ANALYSIS 

Our purpose here is to show that the superstatistical RMT is suitable for the analysis of 
chaotic systems undergoing a process leading to the breaking of a discrete symmetry. In our 
view, symmetry breaking has something in common with the transition from integrability 
to chaos. The presence of a symmetry favors the particular bases in which the eigenvectors 
of the symmetry are components of the eigenstates of the Hamiltonian (e.g., the isospin 
wavefunctions in the nuclear physics problem). In this representation, the Hamiltonian 
matrix is block diagonal and its spectrum consists of independent sequence of eigenvalues. 
As symmetry breaking interactions increase, the eigenstates involved in the different blocks 
mix filling the "empty" places in the Hamiltonian matrix. At a certain stage of the symmetry 
breaking transition, the bases in which the eigenvectors of the symmetry are components 
of the eigenstates of the Hamiltonian loose their special status. The joint matrix-element 
distribution becomes base independent. At this point, we expect superstatistical RMT 
to become suitable for describing the symmetry breaking process with the superstatistical 
parameter v measuring the number of effective matrix elements responsible for symmetry 
breaking. 

In the quartz crystal experiment the authors noted that the block with conserved flip 



2a,l2q|- The initial 



symmetry has much in common with a scalar pseudointegrable system 
state of the transition can be described by an independent superposition of two independent 
semi-Poissonian sequences of equal densities. Applying the method given in Mehta's book 
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ljj , we obtain for NNS distribution 



P2Pi(s) = ^e- 2s (1 + As + 2s 2 ) . (14) 

We have found out that the least square difference between P2pi(s) and the superstatistical 
distribution ffTUj) in the spacing interval < s < 3.0, 

f 30 2 

/ \P2Pl{s) -Pinv X 2 (^ s )\ ( 15 ) 
JO 

has a minimum at v — —0.210. In the following, we show that NNS distributions of reso- 
nances in the quartz crystal experiment [l], the authors noted that the initial state of the 
transition can be described by the distribution ([TO]) allover the symmetry breaking transition 
by allowing v to vary in the range of —0.2 < v < oo. 

In figure 1, we compare the experimental results of the acoustic resonances measured by 
Ellegaard et al. |l] with and the superstatistical NNS distributions corresponding to the 
X 2 , the inverse-^ 2 , an the log- normal distributions of the superstatistical parameter. We 
also show in Fig. 1 the results of a previous analysis of the same data ( figure 5 in Ref. 
js]) with a random matrix model in which assumes that the spectra is composed of three 
independent components, two pseudo-integrable sequences for the conserved symmetry and 
one GOE sequence for the broken symmetry. The best-fit values of the parameters are 
given in Table 1. We quantify the quality of the fits by their absolute average deviation, 
A = jf 2 \Pexp — Pcai.\ i where Pex P , Pcai, and Nl are the experimental, the calculated 
NNS distributions and the number of measures values. The results suggests the validity of 
the superstatistical distribution, even in the initial stages of the breaking of the symmetry. 
The table clearly shown that the NNS distribution obtained from the inverse-^ 2 distribution 
(jlj) agrees with experiment data better than the other distributions even in the initial stages 
of the symmetry breaking process. 
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0.034 


3.387 


0.063 


-0.112 


0.035 


1.100 
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0.037 
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0.061 


-0.025 
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0.047 
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0.046 
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1.064 


0.033 
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0.048 
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0.029 
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0.021 


0.779 


0.025 


1.4 


0.056 


6.67 


0.023 


4.385 


0.024 


0.600 


0.023 


1.7 


0.036 


10.528 


0.020 


7.528 


0.019 


0.469 


0.019 


10 


0.029 


12.183 


0.026 


9.376 


0.026 


0.432 


0.026 



Table 1: The best fit parameters and absolute average deviations A for different radii r 
of the removed corner of the quartz block in the comparison of NNS distributions of acoustic 
resonances and the superstatistical RMT with a y 2 , inverse-^ 2 and log normal parameter 
distributions. The results of analysis of the three- level-sequence model by El-Hady et al. js] 
are also shown 

SUMMARY AND CONCLUSION 

We have described the symmetry breaking of the acoustic resonance in a quartz blocks, 
using a superstatistical model that has been successfully applied to describe systems with 
mixed regular-chaotic dynamics within the framework of RMT. Superstatistical RMT is 
a base independent approach and is suitable to model symmetry breaking only when the 
symmetry representations become mixed enough. Superstatistics arises by allowing the mean 
density of states to fluctuate according to given distribution. We examined three possible 
parameter distributions, namely the x 2 , the inverse-^ 2 and the log normal distributions. Our 
analysis shows that the inverse-^ 2 distribution agrees with experimental spectra of acoustic 
resonances better than the other two, in the same way as in a previous analysis of stochastic 
transition of mixed microwave billiards. The superstatistical parameter v that characterizes 
the fluctuation of the mean level density will hear measure the degree of symmetry breaking. 
We also show that NNS distributions with an inverse-^ 2 superstatistics provide a reasonable 
description of experiment data not only when the system approaches the state of chaos, but 
also in the initial stage of the symmetry breaking transition when base invariance is not 
expected. 
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Figure Caption 

Figure 1. NNS distributions for different radii r of the octant removed from the quartz 
blocks. The experimental data reported in [1] are shown by histgrams. The solid lines 
are the superstatistical results calculted with an inverse-;^ 2 parameter distribution. Results 
for the x 2 an d log normal parameter distributions are shown by dotted and dashed dotted 
dotted crves, respectively. The dashed curves are calculated with the three-level-sequence 
model by El-Hady et al. js]. 
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